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Abstract 

■^ I The Hamiltonian for a particle constrained to move on the surface of a curved nanotube 

^D ' is derived using the methods of differential forms. A two-dimensional Gram-Schmidt 

^:^ ■ orthonormalization procedure is employed to calculate basis functions for determining 

. I the eigenvalues and eigenstates of a tubular arc (a nanotube in the shape of a hyperbolic 

^ I cosine) with several hundred scattering centers. The curvature of the tube is shown to 

'■^ ■ induce bound states that are dependent on the curvature parameters and bend location 

CN ! of the tube. 

T-H . 

PACS number(s): 03.65.Ge, 68.65.-k 
^ ■ Keywords: curved nanotubes, constrained systems 

^: 
!>■ 

O ' 1. Introduction 

O ■ The quantum mechanics of reduced dimensionality systems [1-7] has become a subject of 

Q ■ importance in modelling devices for which the geometric properties of the device become a 

^ ■ factor in influencing its behavior [8-15]. Carbon nanotubes are such objects; the electronic 

'p'' properties of straight nanotubes are determined primarily by the chirality of the tubes and 

p^ . are well understood [16]. However, bent, curved, and toroidal CNT's have been observed 

S ' [17], and their potential application as device elements makes it necessary to model such 

D^' structures as well. 
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In this paper we perform a curved manifold Schrodinger equation calculation to determine 
the spectrum and eigenfunctions for a particle constrained to move on the surface of a 
curved nanotube with delta function site potentials. The methods of differential forms are 
used to derive the Hamiltonian of a particle constrained to move on the surface of a curved 
nanotube. Initially the particle is taken to be in three-dimensional space; it is then confined 
to move on a two-dimensional curved manifold by a potential everywhere normal to the 
surface. This reduction in dimensionality yields a curvature-dependent potential Vd [1- 
7,18-20] that is added to the Hamiltonian. Then, using a hyperbolic cosine shape function, 
the eigenvalues and eigenvectors for the low-lying states of a tubular arc with and without 
(5-function potentials are calculated with a basis set expansion wherein a two-dimensional 
Gram-Schmidt procedure is employed to build the basis states. 

The remainder of this paper is organized as follows: in section 2 the method used to derive 
the Hamiltonian is explained in detail. In section 3, a brief overview of the model used for 



the tubular arc is given. A periodic delta function potential is also introduced to mimic 
atomic sites or defects. In section 4 the method used to solve the Hamiltonian derived in 
section 2 is briefly described. Results for the eigenvectors and eigenvalues of the arc without 
delta functions, as well as the lowest energy eigenvalues with delta sites, are presented. The 
conclusions are given in section 5. 

2. Derivation of the Hamiltonian 

Consider a quantum particle in the neighborhood of a two-dimensional manifold S imbed- 
ded in R^. Any point in the neighborhood of S can be given by the Monge form plus a normal 
term 

r{qi, ^2, qs) = x{qi, gz) + ^363 (1) 

where x describes the manifold, and 63 is the unit vector normal to the surface. Applying 
the exterior derivative operator [21] d gives the one- form 

df=dx + dq^es + qsde^ (2) 

where dq^ is the incremental displacement along the normal to the surface, dx lies in the 
tangent plane and is given by 

dx = (TiCi + 0-262 (3) 

ei and 62 are locally orthonormal unit vectors tangent to the surface; CTi and o"2 are one-forms 
on S. The exterior derivative of any zero- form (scalar function) ip in the neighborhood of S 

is given by 

, > dip . dip . dip , . , , , 

dip = T^dqi + -— dg2 + T^dq^ = ViTi + Z/2T2 + u^dqs (4) 

dqi dq2 Oq^ 

Ui, z/2, and z/3 are zero-form functions and ri, T2, and dq^ are the one-forms in the neighbor- 
hood of the surface S. 

Applying the Hodge * operator to the one-form of Eq. (4) gives 

*dip = iyiT2dq3 + i'2dqzTi + 2/3X1X2 (5) 

A second application of the exterior derivative operator yields 

d* dip = {Aip)TiT2dq3 (6) 

where A is the Laplace-Beltrami operator. The time independent Schrodinger equation 
becomes 



h' 



Aip{r) + V{r)ip{r) = Eip{r) (7) 



2m* 
In the development above, no constraint has been imposed on the particle. To bring 



the particle to the surface, we imagine an oscillator potential V{q'i) = }^m*w'^q^ everywhere 



normal to S acting to enforce the condition that the particle remains in the neighborhood of 
S. In the limit that w -^ oo, then ^3 — * 0. However, it is not enough simply to include this 
term in Eq. (7); to preserve the norm as the particle approaches the surface, it must further 
hold that for gs -^ [2,3,4] 

\'4^{qi,q2,q'i)\'^FdSdq^-^ \x{.qiA2,qz)?dSdq^ (8) 

which allows the identification 

^ = ^ (9) 



where 

F = 1 + 2q^H + qlK 

H and K are the mean and Gauss curvatures respectively and only depend on the surface 
coordinates gi and g2- Performing the standard factorization of x into tangential and normal 
parts and taking the gs — * limit gives two equations when Eq. (9) is substituted into Eq. 
(7), 

--^A.vl/ - -^(ij2 - K)^ = s^ (10) 

2m* 2m* 

and 

The derivative free term in Eq. (10) is the distortion potential Vd, dependent on the mean 
curvature H and the Gauss curvature K of the surface 

VD{qi,q2) = -^{H'-K) (12) 

It should be noted that Vd is not necessarily the only modification to the Laplace-Beltrami 
operator. Even for surfaces possessing symmetry, factors related to the metric can appear in 
the kinetic-energy operator and in general cannot be transformed away [22]. 

To apply the formalism described above to nanotubes, consider a tube of radius a and 
let u be the coordinate axis along i. The surface of the tube can be described by 

x{9, u) = {u — ajS cos 6)i + (/(m) + aa cos 6)j + (a sin 6)k (13) 

where f{u) is the shape function for the axis of the tube, a = (1 + /^)^^^^, /3 = (l + /^)^^/^/„, 
and /„ = dufiu). Applying d to Eq. (13) and using Eq. (3), we find that 

o"i = adO (14a) 

^2 = A(^, s)ds (146) 



and 

Ci = /9(s) sin 6 i — a{s) sin 6 j + cos 6 k (15a) 

e2 = a{s)i + P{s)j (156) 

63 = —P{s) cos9 i + a{s) cos 9 j + sin6' k (15c) 

with the integration measure, axis curvature, and arclength given by A(^, s) = l — aK,{s) cos6, 
K = (1 + fu)~^^'^fuu, and ds = Jl + f^du, respectively. Any point in the neighborhood of 
the tube is given by Eq. (1). The Laplace-Beltrami operator can thus be written 

1 92 91nA (9 1 (9^ (9 In A (9 (9^ dln{Afi) d 

~ ^ W ^ d9 m'^ A2 W ~ ds d^^^d^^ dqs dis ^^^' 

Here 

H = a{l + q^Ki) (17a) 

A = A(l + q3K2) (176) 

where /ti and K2 are the principal curvatures of the tube. The Gaussian curvature K and 
mean curvature H are defined in terms of the principal curvatures by the relationships 

K = K1K2 (18a) 

H = ^1±^ (186) 

For our surface parameterizations the principal curvatures are found to be 

1 , , 

Ki = — (19a) 

a 

1 — a/T:(sj COSC7 
In the limit 53 ^ 0, the kinetic energy operator reduces to 

with the distortion potential given by 

Vnie, s) = -^:r-2u h — W2 (21) 

8m*a"^ [1 — anys) cost^J^ 

Thus, for a quantum particle constrained to move on the surface of a tube, the surface 
Hamiltonian becomes 

ij(^) = -^A, + V^(^,.) (22) 
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and 

is the Hamiltonian due to the normal term. 

3. The Model 

The distorted parts of nanostructures often occur over relatively small sections of the 
object, so as a first step we choose to model a finite tube with pronounced curvature only 
over a small region with hard wall boundary conditions at each end. Further simplification 
follows if we are able to explicitly parameterize by the arclength. As will be seen in the 
next section, this will facilitate performing the large number of integrations that must be 
carried out when computing the Gram-Schmidt coefficients and the surface Hamiltonian 
matrix elements. With this in mind, we choose as our shape function 

/(■u) = cosh kq{u — uq) (24) 

where Kq is the curvature parameter and uq is the turning point of the shape function along 
the u coordinate axis (see Figure 1). The shape function of Eq. (24) has the advantage 
that the arclength as a function of u is one-to-one and analytically invertible. In terms of 
arclength, the curvature of the axis is given by 

<^) = -1 — r-? TT? (25) 

1 + [ko(s - So)]' 

where sq is the turning point along the arclength. Thus kq is the magnitude of the curvature 
of the axis at the turning point. 

As a crude model of atomic sites of the nanotube or as model of defects of a quantum 
waveguide, we add periodic (5-function site potentials to the Hamiltonian 

Na Nr 

V{e, s) = -Ao E E ^(^ - ^jk)S{s - Sk) (26) 

j=i fc=i 

where Aq is the strength of the potential, 6jk is the value of the angle at the j'th site on 
the kth ring, and Sk is the value of the arclength at the A;th ring. A^^ and Nr are the number 
of angular sites per ring and the number of rings, respectively. 

4. Solution Method and Results 

There are many techniques at our disposal for arriving at the eigenvectors and eigenvalues 
of curved nanostructures. In a recent work [22], the eigenstates of a quantum particle con- 
strained to move on the surface of a torus were found using a power series technique. That 



procedure could be adopted in [22] for two reasons, both of which are consequences of a high 
degree of symmetry. First, the distortion potential was only dependent upon one coordinate. 
Secondly, the kinetic energy operator was separable. However, those symmetries are not at 
play here; the distortion potential is a function of both coordinates and the kinetic energy 
operator is not separable. Thus, an alternative method must be employed. The alternative 
chosen here is a basis set expansion. 

The main difficulty faced when trying to expand \E' in a complete set is that, to our 
knowledge, there is no known set of functions {^pk} orthogonal over S, i.e., where 

<Jia2f*{9, s)ipk{0, s) = 6jk (27) 

s 

holds true. To overcome this difficulty, we make use of the Gram-Schmidt procedure [23] 
extended to two dimensions. We choose as our original basis the set of (2M + 1)A^ functions 

im8 , ■ "^^-^ 



U0,s)=Un{O,s) = e'^'sm 



L 



where L is the length of the tube, m = —M, . . . , —1, 0, 1, ... , M , n = 1, . . . ,N, and \j = 1) = 
\m = -M,n = 1), \j = 2) = \m = -M + l,n= 1),..., \j = {2M + 1)N) = \m = M,n = N). 
The orthonormal set {^Pj} is constructed from the non-orthogonal set {C,j}- 

Since we have a legitimate basis, the algorithm is straightforward. The Schrodinger 
equation is solved as a matrix eigenvector-eigenvalue equation with the eigenvectors being 
the coefficients of the expansion and the elements of the Hamiltonian matrix are defined by 

H^;^ = f a,a,^*{9, s)H^'^^MO, s) (28) 



s 

Since the Hamiltonian in Eq. (22) is invariant under 6 ^ —6 we expect the solutions to 
segregate themselves into positive and negative parity solutions. This will prove useful as a 
check on the reliability of the numerical results. It should also be noted that in the presence 
of delta functions, we expect the number of angular peaks in the ground state probability 
density to be equal to the number of delta sites per ring. This will also serve as a useful 
check on the numerics. 

Using m* = rrie the e < results for the tubular arc without delta function site potentials 
are presented in tables 1-4. As can be seen, all e < 0, i.e. m = states, are of positive 
parity. The |m| > states alternate between positive and negative parity. Table 1 shows the 
results for the ground state of tubes with several curvature parameters and turning points. 
Convergence was achieved with 5-digit accuracy using a 20-state (M = 2, A^ = 4) basis set 
expansion. When k,q = 0, the tube is straight and the system can be solved analytically; the 
distortion potential is a constant that can be subtracted from the eigenvalues but is left here 



for comparison with the arc. It can be seen that the energy of the ground state is lowered 
when the value of the curvature parameter is increased. However, curvature has little effect 
on the energies of the excited states (tables 2-4). It should also be noted that the energy of 
the ground and excited states are sensitive to the position of the bend as well. If the turning 
point is closer to either end of the tube than it is to L/2, the eigenenergies of the ground, 
first and third excited states are again lowered, however, the energy of the second excited 
state is raised. 

Curvature and bend locations have an effect on the charge (probability) density as well. 
Although the charge density of the straight tube is angular independent, curvature induces 
angular dependence in the charge density of the ground state that is maximum at the point 
of minimum radius of curvature {6 = tc,s = sq), i.e. the point of maximum curvature. The 
charge density of the ground state at ^ = tt is shown in Figure 2. There is no angular 
dependence in the excited states, but curvature and bend location still have an effect on the 
charge density. The nodes and peaks of the density are shifted and the height of the peaks 
are no longer even (Figs. 3-5). 

The eigenvalues of the arc with 1170 {Na = 6,Nr = 195) delta function potentials 
arranged in an armchair configuration are given in table 5. These results were obtained 
with a 52-state (M = 6, A^ = 4) expansion. The value for the strength of the potential was 
chosen to be of the same order of magnitude as the deepest part of the distortion potential 
well. As can be seen, the ground state is more sensitive to the curvature than the excited 
states. The slight increase in energy for the third excited states shown in the table may be 
due to a need for more basis states in the expansion. In the presence of the delta function 
potentials, charge localization still occurs in the region of curvature (Figure 6); however, 
there are several islands of localization, with the maximum occurring near the center of the 
potential well. As expected, the number of islands is equal to the number of delta function 
sites in each ring, i.e., if there are 6 deltas per ring, there are 6 peaks. 

5. Conclusion 

The eigenstates and eigenvalues of a particle constrained to move on the surface of a 
tubular arc with and without delta function sites were computed for several values of k,q 
and Sq. It was shown that the energy and density were not only sensitive to the strength 
of the curvature but also to the location of the bend. There was strong interplay between 
the curvature and delta function potentials demonstrated by the density plot of Figure 6. 
Charge gets localized by the deltas as well as curvature. The density is still peaked over 
s = Sq, but the maximum peak is no longer at ^ = vr, the center of the well. This peak is 
slightly displaced due the deltas competing with the curvature as well as each other. What 
this type of interplay would mean for electron transport is an open question and a topic 
currently under investigation. 
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Table 1: Ground state eigenfunctions and eigenvalues for the tubular arc without delta 
function site potentials. For all tubes considered here a = O.Sbnm, L = lOOnm. Coefficients 
not listed are at least an order of magnitude smaller than those given. 



Ko 



So 



"^mn, d = 0.85nm, L 



lOOnm 



e{meV) 



0.00 
0.75 

0.75 

0.95 

0.95 

1.00 
1.00 



51.87 
55.60 
52.37 
57.08 

52.50 

57.45 



1.15 73.79 



^ 



.0612 sin 



01 = .uuizbmi^Y^) 
^01 = -.0022 cos ^sin(f^ 

^^'^ - .0253sin('IS^ 
01 = .0023 cos ^ sin (f^ 
.0019 cos ^sinff^ 



+.0021 cos ^sin(f§) 



+ .0555 sin 

37rs 
100 



.0042 Sim 



^01 



ns 
100 

.0040sin(f 
.0541sin(f^ 

37rs 



' TTS ' 

50' 



.0135sin(f) 



100' 



.0227sin(ff ) + .0014 cos ^sin(f) 



-.0024 cos ^ sin (f§) 



.0026 cos^ sin(f^) - .0526 sin(f^) + .0061 sin(f ) 
+ .0302sin(f§) - .0066sin(f ) 
■01 = -.0027 cos ^sin(f^) + .0500sin(f^) + .0012 cos ^ sin 
-.0179 sin(f ) + .0021 cos^ sin(f^) - .0248 sin(ff^, 
+.0183sin(f) 

*oi = .0027 cos ^sin(f^) - .0517sin(f^) + .0065sin(f ) 
-.0025 cos ^ sin (Iff) + .0315sin(f^) - .0074sin(f ) 
^01 = -.0028 cos ^sin(^) + .0488sin(^) + .0013 cos ^ sin 



.0117sin(ff) 



TTS ' 

50' 



-. 0020 COS ^sin(f; 



TTS ' 

50' 



.0250sin(^) - .0022 cos ^ sin f^ 



100 



100' 

-.0189sin(f ) + .0021cos^sin(f^) 

+.0202sin(f ) 

^02 = .0027 cos ^sin(f^) - .0356sin(f^) - .0036 cos ^sin(f) 

+.0438sin(f ) + .0022 cos ^ sin (f^) - .0235sin(ff ) - .0048sin(f) 



25' 



-13.1423 
-13.4068 

-13.4262 

-13.5793 

-13.6313 

-13.6522 

-13.7191 

-14.2313 



Table 2: Eigenfunctions and eigenvalues of the first excited states of the tubular arc. Coef- 
ficients not listed are at least an order of magnitude smaller than those given. 



0.85nm, L = lOOnm 



Kq 



So 



^ mm "■ 



e(meV) 



0.00 
0.75 
0.75 
0.95 
0.95 
1.00 
1.00 
1.15 



51.87 
55.60 
52.37 
57.08 
52.50 
57.45 
73.79 



^02 
^02 
^02 
^02 
^02 
^02 
^02 
^01 



.0612sin(f) 
.0067sin(fii; 
.0185sin(f^; 
.0103sin(f^; 
.0257 sin(^; 

TTS 

aoo 
-.0275sin(fi^: 

TTS 



100 

-.0114sin(^ 



+ .0606sin( 
+ .0577sin( 
+ .0600sin( 
+ .0545 sin( 



■KS_' 

50 
50 

■KS 

50 

■KS ' 

50 



.0488sin(f^) 



-.0597sin(f 
-.0536sin(f 



+ .0045sin(|f) 
+ .0084sin(ff) 
+ .0061sin(ff) 
+ .0098sin(f§) 



-.0025sin( 



25- 



.0067sin(ff) 
.0100sin(ff; 



-.0038sin(f ; 



.0011 sin 

— nso/ •-^---^^100 7 i .0041 sin, 
.0365sin(f ) - .0057sin(f^) - .0012sin(f 



TTS ' 

25' 

TTS ^ 

25^ 



-13.0296 

-13.0308 

-13.043 

-13.0308 

-13.0476 

-13.0312 

-13.0497 

-13.1016 
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Table 3: Eigenf unctions and eigenvalues of the second excited states of the tubular arc. 
Coefficients not listed are at least an order of magnitude smaller than those given. 



0.85nr7i, L = lOOnm 



Kq 



So 



^ 



.0612sin(f§ 



mm 



a 



e{meV) 



0.00 
0.75 
0.75 
0.95 
0.95 
1.00 
1.00 
1.15 



51.87 
55.60 
52.37 
57.08 
52.50 
57.45 
73.79 



^03 
^03 
^03 
^03 
^03 
^03 
^03 
^03 



,0249sin(f^) - .0068sin(f ) + .0554sin(f^) - 
-.0211 sin(f^) + .0151 sin(f ) - .0546 sin(f^ 
-.0295sin(f^) + .0103sin(f ) - .0524sin(f^ 
-.0225sin(f^) + .0208sin(f ) - .0516sin(f^ 
.0306sin(f^) - .0115sin(f ) + .0514sin(ff ) - 
-.0227sin(f^) + .0221sin(f ) - .0508sin(f^ 
-.OlOOsin(f^) + .0222sin(f ) + .0559sin(f^ 



-.0035sin(f ) 
+ .0093sin(f ) 
+ .0049sin(f ) 
+ .0012sin(f ) 

-.0053sin(f ) 
+ .0013sin(f ) 
+ .0041sin(f) 



-12.8416 
-12.9331 
-12.9169 
-12.9478 
-12.9182 
-12.9515 
-12.9178 
-12.874 



Table 4: Eigenf unctions and eigenvalues of the third excited states of the tubular arc. Coef- 
ficients not listed are at least an order of magnitude smaller than those given. 



Kq 



So 



\1' 

^ r 



0.85nm, L = lOOnm 



e{meV) 



0.00 
0.75 
0.75 
0.95 
0.95 
1.00 
1.00 
1.15 



51.87 
55.60 
52.37 
57.08 
52.50 
57.45 
73.79 



^04 
^04 
^04 
^04 
^04 
^04 
^04 
^04 



.0612sin(f) 
-.0022sin(f^) 



.0048sin(§f) 



.0065 sim 
.0034 Sim 
.0093 sim 
.0038 sim 
.0102 sim 
.0012 sim 



• TTS ' 

100' 

• ns ' 

100' 

• ns ' 

100' 

• ns ' 

100' 

• ns ' 

100' 

• ns ^ 

100' 



-.0027 sin 

- .0074 sin 

- .0049 sin 
+ .0013 sim 

- .0056 sin 

- .0028 sin 



3ns - 

100 

' ns\ 

50'' 

50 ' 
ns\ 

50'' 
ns\ 

■50'' 



+ .0610sin(f ^ 
.0133sin(f^)- 
-.0606sin(f ) 



-.0194 sin 
-.0083 sim 
- .0211 sin 
+ .0057 sim 



' 3ns ^ 

100' 

3ns ^ 

100^ 

• 3ns ^ 

100' 
3ns ^ 

100^ 



.0593sin(f ) 



.0571 
.0605 
.0562 
.0609 



sin(f) 
sin(f) 
sin(f) 
sin(f) 



-12.5784 
-12.5799 
-12.6136 
-12.5805 
-12.6398 
-12.5804 
-12.6475 
-12.25592 



Table 5: Eigenvalues of the tubular arc with 1170 {Na = 6,Nr = 195) 5-function site 
potentials arranged in an armchair configuration with Aq = 400meV ■ nm. The subscripts 
on the eigenvalues are not quantum numbers. They refer to the order of the energies: is 
the ground state; 1 is the first excited state and so on. 



Ko Sq 


€o{meV) 


ei{meV) 


e2imeV) 


esimeV) 


0.00 — 
0.95 52.37 
1.00 52.50 


-872.208 

-872.587 
-872.655 


-872.095 
-872.115 
-872.119 


-871.907 
-872.012 
-872.018 


-871.643 
-871.627 
-871.624 
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Figure 1: A tubular arc. 
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Figure 2: Probability densities of the states given in Table I at ^ = vr. The color code is 
as follows and is the same for the remaining figures: solid black (kq = 0.00); dashed blue 
(kq = 0.75, So = 51.87); dashed turquoise (kq = 0.75, Sq = 55.60); solid turquoise (kq = 
0.95, So = 52.37); solid blue (kq = 0.95, Sq = 57.08); dashed pink (kq = 1.00, Sq = 52.50); 
dashed green (kq = 1.00, sq = 57.45); solid pink (kq = 1.15, sq = 73.79). 
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Figure 3: Probability densities of the states given in Table II at = tt. 
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Figure 4: Probability densities of the states given in Table III at ^ = tt. 
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Figure 5: Probability densities of the states given in Table IV at ^ = tt. 
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Figure 6: Probability density of the ground state of a tubular arc (kq 
with delta function site potentials. 



1.00, So = 57.45) 
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